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ABSTRACT 

An aeroacoustic model test has been conducted to 
investigate the mechanisms of sound generation on 
high-lift wing configurations. This paper presents an 
analysis of flap side-edge noise, which is often the most 
dominant source. A model of a main element wing 
section with a half-span flap was tested at low speeds of 
up to a Mach number of 0.17, corresponding to a wing 
chord Reynolds number of approximately 1.7 million. 
Results are presented for flat (or blunt), flanged, and 
round flap-edge geometries, with and without 
boundary- layer tripping, deployed at both moderate and 
high flap angles. The acoustic database is obtained 
from a Small Aperture Directional Array (SADA) of 
microphones, which was constructed to electronically 
steer to different regions of the model and to obtain far- 
field noise spectra and directivity from these regions. 
The basic flap-edge aerodynamics is established by 
static surface pressure data, as well as by 
Computational Fluid Dynamics (CFD) calculations and 
simplified edge flow analyses. Distributions of 
unsteady pressure sensors over the flap allow the noise 
source regions to be defined and quantified via cross- 
spectral diagnostics using the SADA output. It is found 
that shear layer instability and related pressure scatter is 
the primary noise mechanism. For the flat edge flap, 
two noise prediction methods based on unsteady- 
surface-pressure measurements are evaluated and 
compared to measured noise. One is a new causality 
spectral approach developed here. The other is a new 
application of an edge-noise scatter prediction method. 
The good comparisons for both approaches suggest that 
much of the physics is captured by the prediction 
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models. Areas of disagreement appear to reveal when 
the assumed edge noise mechanism does not fully 
define the noise production. For the different edge 
conditions, extensive spectra and directivity are 
presented. Significantly, for each edge configuration, 
the spectra for different flow speeds, flap angles, and 
surface roughness were successfully scaled by utilizing 
aerodynamic performance and boundary layer scaling 
methods developed herein. 

SYMBOLS 

a 0 medium speed of sound 
c flap chordlength 

C N normal force coefficient with respect to c 

C static pressure coefficient 

COP s coherent output power spectrum of unsteady 
surface pressure with respect to far-field noise 
d distance from one sensor to another 
D directivity factor, Eq. (13) 

dS{ y) elemental surface area at y 

/ frequency 

/ 1/3 one-third octave band center frequency 

A f spectrum frequency bandwidth 

G i{ auto-spectrum of noise measured by SADA 
G s auto-spectrum of unsteady surface pressure at 

sensor 

G a y cross-spectrum between outputs of SADA and 

surface pressure sensor 
i pressure sensor location number 

j V-T 

k acoustic wave number = col a 0 

correlation length scale in chordwise edge 
direction 

^3 correlation length scale in spanwise direction 
from edge 

L length of chordwise section that a sensor 

represents 

L' lift per unit span 

M v convective Mach number, U c / a Q 
M c average A/ , see Eq. (8). 

A/ 0 tunnel Mach number. U 0 / a () 

n normal vector to surface at y 
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R c Reynolds number based on c and U 0 

p a acoustic pressure time history 

p s surface pressure time history 

P a Fourier transform of p a 

P s Fourier transform of p s 

q c dynamic pressure based on convective speed 

Uc 

q 0 tunnel dynamic pressure 

r distance = |r| 

r effective source-to-observer distance for 

quiescent field radiation, see Fig. 19 
r vector distance = x - y 
SADA small aperture directional array 
t time 

u su inboard velocity along suction side 

u pr inboard velocity along pressure side 

u v velocity at radius r 0 from center of vortex, see 
Fig. 9(c) 

U c convection velocity 

U h hydrodynamic convection speed from the 

pressure to suction side, see Fig. 13 
U' h hydrodynamic convection speed over the 

suction side, see Fig. 13 

U flow velocity at edge on pressure side, Eq. (4) 

U s flow velocity at edge on suction side, Eq. (4) 

U 0 tunnel velocity 

x chordwise distance from the flap leading edge 
x noise observer location vector 

x' effective observer location vector for radiation 
in quiescent field, see Fig. (19) 
y spanwise distance from the flap edge 

y surface noise source location vector 

z height above surface sensor, see Fig. 8 
a flap angle with respect to the main element 
/3 C edge convective-flow skew angle, see Fig. 8 
y circulation density 

yl s coherence function, see Eq. (9) 

r vortex circulation 

8 boundary layer thickness 

<5 0 boundary layer thickness at airfoil zero angle 

of attack 

g coherence decay factor for £ 3 

q coherence decay factor for 

9 angle between n and x, see Fig. 19 

9’ angle between n and x', see Fig. 19 

& observer azimuth angle defined for Eq. (2 1 ) 
p medium density 

r retarded time, Eq. (12) 

r a s noise transmission time from sensor to SADA 

<p SADA elevation (flyover) angle, see Figs. 5 
and 24 


(p as cross-spectral phase between SADA and 
sensor outputs 

t | f SADA azimuth angle, see Fig. 24 

ft) radian frequency = 2 7tf 

INTRODUCTION 

Airframe noise can be dominant during airport 
approach and landing when the engines are at low 
power and the high-lift systems and landing gears are 
deployed. This becomes particularly true as present- 
day propulsive systems become quieter 1 . As a result, 
there has been an increased emphasis placed on the 
measurement and modeling of non-propulsive 
components such as flaps, slats, and undercarriage. 

2 

As reviewed by Crighton", a number of studies of 
airframe noise were conducted in the 1970's and early 
1980's. An early evaluation was performed by Hardin . 
Empirical airframe noise studies and prediction 
developments include those of Fink 4 and Fink and 
Schlinker . A series of airfoil self-noise experiments 
were performed by Brooks and Hodgson 6 and Brooks 
and Marcolini * ’ for trailing edge noise and wing tip 
noise. The results of these studies formed the basis of a 
comprehensive self-noise prediction method 10 for 
isolated airfoils. As part of a wing and flap high-lift 
system, the flap is much more loaded aerodynamically 
than it would be if isolated. Because of this, it has been 
found capable of producing much more intense noise. 
Block 11 in wing, flap, and landing gear interaction 
studies found flaps to contribute significantly to the 
overall noise. Kendall 12 and Kendall and Ahtye 13 , 
using an elliptical acoustic mirror, found strong 
localized flap edge noise. This was confirmed by Fink 
and SchlinkeE in component interaction studies. 
Mclnerny et al. 14 , Ahtye et al. 15 , and Miller and 
Meecham performed cross-correlation studies 
between unsteady surface pressures and noise field 
measurements for the tip region of an isolated wing, 
single slotted flap, and triple slotted flaps, respectively. 
The side edges of the multiple flaps were found to 
significantly exceed other airframe noise sources 16 . 

The 1990's produced an increase in airframe noise 
research activity 17 , particularly due to the NASA 
Advanced Subsonic Technology (AST) program. 
Several tests are particularly notable. A 4.7% scale 
DC- 10 aircraft model was tested in the NASA Ames 40 
by 80 foot wind tunnel, as reported by Bent et al. 18 , 
Hayes et al. and Guo et al. ' . Inflow microphones, 
a phased-microphone array, and a parabolic mirror 
directional microphone system were used along with 
unsteady surface pressure sensors on inboard and 
outboard flaps. The flap edge noise was found to 
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dominate other noise sources. Significant correlations 
were found between edge pressures and the measured 
noise . Noise reduction concepts were evaluated - . A 
series of tests of a large unswept wing (2.5 ft. chord) 
and half-span Fowler flap were conducted in the NASA 
Ames 7x10 foot wind tunnel, as reported by Storms et 
al. 23 , Horne et al. 24 , and Storms et al. 25 . The tests 
provided basic aerodynamic data and, although the 
tunnel was hard-walled, limited acoustics were obtained 
using large phased arrays of microphones. A 
computational study by Khorrami et al. 26 provided 
substantial agreement with the data. This was used to 
examine two possible noise source models, namely, a 
vortex-instability model and a shear layer vortex-sheet 
model. 

The present paper concerns a wing and flap model 
tested in the Quiet Flow Facility (QFF) at NASA 
Langley. The model is a NACA 63 2 -215 wing with a 
30% chord half-span Fowler flap. This is the same as 
that used in the aforementioned 7x10 foot wind tunnel 
test at NASA Ames, except here the model is about one 
half the size. As reported by Macaraeg 17 , this model in 
the QFF has provided the means to closely examine the 
aerodynamic and acoustic physics for slats and flaps. 
Measurements of the flow field in the QFF, by 
Radezrsky et al. 27 included hot-wire, hot film, 5-hole 
probe surveys, laser light sheet, and flap surface oil 
flows. These measurements revealed a dominant flap 
vortex structure resulting from the merger of two 
upstream vortices - one strong vortex, formed from the 
pressure side to around the flap edge, and a weaker 
vortex formed at the flap side edge on the suction side. 

In the vicinity of the trailing edge, the vortex is far 
removed from the flap surface. Computational efforts 

tb* 29 

by Khorrami et al." and Takallu and Laflin" using 
Reynolds Averaged Navier-Stokes solutions (RANS) 
duplicated the key mean features of the edge flow. 
Streett 30 developed a computation framework for the 
simulation of the fluctuating flowfield associated with 
this complex flap-edge vortex system. Streett’s 
computations, utilizing a calculated mean flow field 28 , 
further crystallized the shear layer instability and 
vortex-instability disturbance models" 6 for noise 
production. Linear stability analysis determined 
dominant frequency ranges of unstable flow 
disturbances 31 . Guo 32 33 , in a similar time frame, 
followed with a semi-analytical and semi-empirical 
prediction model of this shear layer instability 
mechanism. Predictions from this model compared 
well with flap edge noise data when certain scale 
parameters were used. 

The initial aeroacoustic measurements for an 
instrumented version of the above model tested in the 
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QFF were presented by Meadows et al. . 
Measurements included flap-edge noise-source location 
mapping by a large directional (phased) microphone 
array system, flap-edge noise spectra and directivity by 
a smaller array, and cross-spectra between unsteady 
surface pressure sensors about the flap edge. Details of 
the microphone array design and methodology used in 
the testing was presented by Humphreys et al. . 
Microphone array testing methodology was refined and 
quantified using the QFF systems, as reported by 
Brooks et al. 36 . The present study builds upon this 
work. 

In this study, the generation and radiation of Hap 
edge noise for the flat (or blunt), flanged, and round 
flap edge configurations are examined. The basic flow 
pattern about the edge is studied using Computational 
Fluid Dynamic (CFD) calculations and measured static 
pressure distributions. Simplified flow calculations are 
then developed to provide key aerodynamic parameters 
needed for noise prediction and scaling. Cross-spectral 
amplitude and phase between unsteady surface pressure 
sensors over the flap edge surface are analyzed to 
reveal the character of the hydrodynamic pressure field 
due to turbulent flow and the near-field flap-edge noise 
generation. Coherent Output Power (COP) spectra 
diagnostics using the measured pressures and the noise 
provide a measure of the noise source distribution along 
the flap edge. The noise source thus determined is 
examined for consistency with the previously 
mentioned shear layer instability mechanism -6 ' 3033 . For 
the flat edge flap, separate noise prediction methods are 
developed and validated from (1) a causality approach 
that connects the noise to the cross-spectra between the 
surface pressure and far- f i eld noise through 
fundamental aeroacoustic formulations and (2) an edge- 
noise scatter solution. Both methods utilize the surface 
pressure measurements on the suction and pressure 
sides near the flap edge. Next the noise spectra and 
directivity are presented for three edge configurations 
for different surface roughness, flap angles, and flow' 
speeds. The spectra are then examined for scalability 
for each configuration using flap mean lift and 
boundary layer thickness descriptions. 

TEST SETUP AND METHOD 

The test model apparatus is shown mounted in the 
Quiet Flow Facility (QFF) in Fig. 1 . The QFF is a quiet 
open-jet facility designed for anechoic acoustic testing. 
For the present airframe model testing, a 2 by 3 foot 
rectangular open-jet nozzle is employed. The model is 
a NACA 63 2 -215 main-element airfoil (16 inch chord 
and 36 inch span) with an attached half-span Fowler 
flap (4.8 inch chord). The flap is attached by an 
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adjustable set of "U" brackets to minimize bracket 
interference with the ideal flap flow field. The model is 
held in place by vertical side plates, which are 
themselves rigidly mounted to side plate supports of the 
nozzle. In the photo, the model is visible through the 
Plexiglas windows located on the side plates. The main 
element airfoil and flap are instrumented with static 
pressure ports and unsteady pressure sensors 34 . 

A view of the main element and flap in the vicinity 
of the flap edge is sketched in Fig. 2. The flat edge flap 
is shown accompanied by edge modifications. When 
attached, the flange edge produces a cavity depth of 1/8 
in. The flange thickness is 0.05 in. The round edge 
attachment is a half-circle cross-section shape that 
matches the airfoil contour. The effect of surface 
roughness on the flap edge noise was examined by 
applying grit. For the flat edges, #60 grit at a density of 
about 70 particles per square inch was applied on the 
edge and both suction and pressure side surfaces over a 
2 in. span. For the round edge, #120 grit at about 800 
particles per square inch was applied, but was restricted 
to one half of the round edge surface area - towards the 
flap's pressure side. The intent of the grit was to 
produce thickened and well-developed turbulent 
boundary layers in the vicinity of the side edge. For 
this paper, the main element angle was set at 16° and 
two flap angles, a =29° and 39°, were tested. The gap 
and overlap settings for these angles are shown in Fig. 
3. The positions of the flush-mounted unsteady 
pressure sensors in the flap edge vicinity are shown in 
Fig. 4. The chordwise distance from the leading edge 
is x and the span wise distance from the side edge is y . 



FIGURE 1. Test apparatus with SADA mounted on pivotal 
boom in QFF. 



FIGURE 2. Sketch of flap edge treatments. 


MAIN 

ELEMENT OVERLAP 



FIGURE 3. Flap and gap geometry. 

As will be discussed, the sensors of present interest are 
those on the pressure and suction surface. These 
sensors are Kulite model LQ-34-064-5A. They are 
aligned spanwise at .06, .81, and 1.81 in. (sections A, B, 
and C, respectively) from the edge. The chordwise 
position for each sensor is given in Table 1 . 

The far-field acoustics of the model are measured 
by Small Aperture Directional Array (SADA), which is 
seen in Fig. 1 to be mounted on a pivotal boom 
positioned by rotational stepping motors. The SADA is 
always 5 ft. from the center of the main element trailing 
edge. It consists of 33 B&K 1/8-inch microphones 
projecting from an acoustically treated metal frame. 
The aperture of the array is small, with a maximum 
diagonal aperture of 7.76 inches. The small size 
reduces bias error by locating all the microphones in the 
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array within approximately the same source directivity, 
regardless of SADA’s elevation or azimuth position 
about the model. In Fig. 5, the SADA measurement 
positions are drawn in a side view (opposite side to that 
of Fig. 1 ) of the test setup. The SADA is shown located 
in a plane perpendicular to and centered on the span of 
the model, corresponding to zero azimuthal angle (\| f = 
0°). The position of SADA in the photo of Fig. 1 
corresponds to an elevation angle 0 = -124° in the 
drawing. In Fig. 5, the SADA is seen positioned at <j) = 
107°, on the pressure side of the model. The open jet 
shear layer boundaries (defined at 10% and 90 % of the 
potential core velocity) are shown as measured along 
the \| f =0° plane. A mean shear line is shown, which is 
part of a curved three-dimensional mean shear surface 
defined mathematically from the shear layer 
measurements. This is used in SADA processing to 
determine shear layer refraction corrections. The 
drawing illustrates the refracted noise ray path from the 
flap edge source region to the microphone. 


<b = -39° 
□ 


-56° 

□ 



□ 

•107° 


□ 

-124° 



FIGURE 5. Sketch of test setup. The noise ray path from the 
flap edge to the SADA is illustrated. 


SUCTION EDGE PRESSURE 



C B A ABC 


FIGURE 4. Drawing of unsteady surface pressure sensor 
distribution. 


Sensor Coordinates 

" 

# 

x(inch) 

# 

x(inch) 

# 

x(inch) 

# 

x(inch) 

10 

0.12 

20 

0.95 

30 

0.12 

39 

3.02 

ii 

0.54 

21 

2.62 

31 

0.95 

40 

2.18 

12 

0.95 

22 

3.45 

32 

1.78 

41 

1.35 

13 

1.37 

23 

3.99 

33 

2.20 

42 

0.93 

14 

1.78 

24 

0.12 

34 

2.62 

43 

4.80 

15 

2.62 

25 

1.78 

35 

303 

44 

4.68 

16 

3.03 

26 

3.45 

i 36 

3 36 

45 

3.02 

, 17 

3.36 

27 

3.99 

1 37 

0.00 

46 

1.35 

1 18 

3.99 

29 

0.00 

; 38 

0,12 

47 

0.93 

1 19 

0 12 

- 


i - 





y A = 0.06 inch y B = 0.81 inch y c = 1 .81 inch 


TABLE 1. Pressure sensor coordinates. 


Data acquisition and post-processing 

The array microphones and surface pressure 
sensors employed acquisition hardware consisting of 
transient data recorders controlled by a workstation. All 
35 microphone channels (including 2 reference 
microhones) were recorded with a 14-bit dynamic 
range, simultaneously with 32 pressure sensor channels 
using a 12-bit range, at a sampling rate of 142.857 kHz. 
Two million 2-byte samples were taken for each 
acquisition. The microphone signals were high pass 
filtered at 300 Hz. All channels had anti-aliasing filters 
set at 50 kHz, which is substantially below the 71.43 
kHz Nyquist frequency. 

Microphone and pressure sensor calibration data 
were accounted for in the post-processing. For the 
SADA microphones, regular pistonphone and injection 
calibrations of amplitude and phase were made. 
Amplitude and phase calibrations for the pressure 
sensors employed a miniature speaker-driver capable of 
high frequency output. The measured outputs were 
referenced to the output of a 1/8 in. B&K Model 4133 
microphone. (The high frequency outputs of the 
present Kulite sensors are unfortunately limited. In this 
report, surface pressure spectral data is limited 
generally to 13.5 kHz, where flat frequency response 
and signal-to-noise are good.) Initial post processing of 
the test data begins with the computation of the cross- 
spectral matrix for each data set. The computation ot 
the individual matrix elements is performed using Fast 
Fourier Transforms (FFT) of the original data 


5 


American Institute of Aeronautics and Astronautics 



ensemble. All data are segmented into 1000 non- 
overlapping blocks each containing 2 12 samples, 
yielding a frequency resolution of 34.88 Hz. A 
Hamming window is used. 

A conventional beamforming approach, employing 
matrices of cross-spectra between the array 
microphones 35,36 , is used to electronically “steer” the 
array to chosen noise source locations. The processing 
accounts for mean amplitude and phase changes due to 
refracted sound transmission through the shear layer to 
the individual microphones of the arrays. A mean 
refracted ray path is illustrated in Fig. 5. The correction 
terms are calculated 33 using Snell’s law in Amiet’s 
method , modified to account for a curved three- 
dimensional mean shear surface defined in the shear 
layer. A key feature of the array processing is that 
spatial resolution (or sensing area over noise source 
regions) can be controlled independently of frequency 
and steering-direction over broad frequency ranges. 
The microphone shading algorithms methodology used 
is adapted from Refs. 38 and 39 and evaluated with 
respect to the present test in Refs. 35 and 36. Note that 
for each test case, the cross-spectral matrix has a 
corresponding background matrix subtracted from it to 
remove extraneous system noise (measured microphone 
and sensor noise for zero tunnel flow speed). The array 
processing references levels to an equivalent single 
microphone measurement. Spectra data are determined 
for each narrowband frequency (34.88 Hz resolution 
bandwidth) of interest. Other spectral bandwidths that 
are presented in this paper are formed from the 
narrowband spectra. 

FLAP EDGE FLOW FIELD 

In this section, the basic flap edge flow is 
examined with respect to parameters required to 
evaluate the unsteady surface pressures and related 
noise field. 

Basic aerodynamics 

Extensive aerodynamic measurements for the 
present model have been reported by Radezrsky et al. 27 . 
The model was shown to function as a high-lift device, 
with the main element and flap properly interacting 
aerodynamically. The elements are close enough that 
the flow acceleration about the leading edge of the flap 
significantly reduces the required pressure recovery at 
the main element trailing edge, but the elements are 
separated sufficiently so that the viscous boundary 
layers do not merge. This increases the overall lift, 
especially on the main element, compared to lifts 
obtainable separately. The QFF facility produces a 


maximum Mach number of 0.17 for this model 
configuration, which corresponds to a main element 
chord Reynolds number of 1.7 x 10 6 , In order to 
maintain attached flow on the flap, the boundary layer 
transition was fixed by serrated tape applied to the 
lower surface of the main element at 30% chord and on 
the leading edge of the flap. Pressure coefficient plots 
revealed very similar performance to the somewhat 
larger Reynolds number conditions of the similar 
model 23 tested in the Ames closed wall 7 x 10 foot 
tunnel. In the QFF, the flap angle with respect to the 
main element was a= 29° and 39°, whereas the main 
element was set at 16° and 20° angle of attack to the 
tunnel centerline. (Note that 16°, for the main element, 
is approximately equivalent to an angle of attack of 
about 5° in the closed wall tunnel.) The flap flow field 
was found to be dictated almost entirely by the flap 
angle, which is measured with respect to the main 
element, and not the main element angle. 

For the present QFF testing, pressure and lift 
distributions for the flap are presented in Fig. 6. The 
main element angle was 16°. The gap and overlap 
settings, shown in Fig. 3, differ only slightly from those 
of Ref. 27. Static pressure coefficient distributions at 
three spanwise locations of the flap are shown in Fig. 6 
for the tunnel Mach number M 0 = 0.17 for the two a 
values. The spanwise cuts are shown for yfc - 0.027, 
0.208, and 1.875. The ratio yfc is the distance from 
the flap edge compared to the flap chord length c. At 
yfc - 1.875, at the center of the flap section, the 
expected two-dimensional lift distribution behavior 
with high suction peaks is shown for both angles. As 
yfc decreases (meaning the flap side is approached), 
the high suction peak at the forward (leading edge) 
stations are reduced and the pressure differential 
diminishes. Near the side edge, a low-pressure region 
exists at a downstream section of the chord, which is 
due to a strong vortex being formed on the suction side. 
Also shown in Fig. 6 is the normal force (normal to 
chordline) coefficient C N , with respect to c, versus 
yfc. An additional yfc location of 0.625 is 
represented here. It is seen that the sectional lift is 
diminished as the side edge is approached except for an 
increase very near the edge due to the presence of the 
strong vortex on the suction surface. At the inboard 
station yfc - 1.875, C N = 1.213 and 1.567 for a = 
29° and 39°, respectively. The ratios of C N and a 
values show almost a linear dependence of lift to flap 
angle. 

The vortex found on the suction surface near the 
flap edge was shown in Ref. 27 to be a result of the 
strong primary vortex and a weaker vortex merging. 
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FIGURE 6. Pressure coefficient distributions and normal force 
coefficient distribution for two flap angles. 

The primary vortex is formed along the pressure side 
(bottom) edge and grows in size in the streamwise 
direction, and a weaker vortex is formed near the 
suction surface edge. Steady RANS computations of 
Ref. 28 found agreement with the basic measured 
features of the merger of the dual vortex system and the 
general location of the resultant vortex. For both the 
experiment and calculations, the vortex bursts above the 
suction side surface for the 39° flap angle case. This 
bursting occurs when the local flow angularity is too 
high or the axial velocity component is too low. Figure 
7 shows portions of the RANS solutions for the two 
flap angle QFF test cases of the present study. The 
contours show lines of constant static pressure on the 
surface. Intervals between the lines correspond to 
intervals in C p of .346. The two component vectors 

shown are the calculated velocities over a projected 
surface defined at 0.035 in. (approximately a boundary 
layer thickness) above the suction and pressure 
surfaces. Only the edge velocity vectors from the 
pressure side are seen because of the oblique view of 
Fig. 7. The flow about the side edge surface is omitted 
for clarity. The vector pattern clearly shows the 
presence of the resultant vortex and its strong influence 
on the flap edge flow field. The vortex is trailed 
downstream of the model, but the vectors show the 
formation of the vortex is essentially attached at the top 
(suction) edge surface. The attachment is seen to be 
just aft of mid-chord for the 29° flap angle case, but 
slightly forward of mid-chord for the 39° flap angle. 
The vortex strength is mostly defined by the strong 


sheared-flow velocity across the pressure surface edge 
which wraps around the vortex and "feeds" it. 


M, 




Of primary interest for this study are flow 
parameters that provide guidance in determining noise 
sources and provide pertinent input to prediction theory. 
If the flap edge noise problem is indeed an edge 
scattering problem, one would view the boundary layer 
character and associated velocities as primary 
parameters. One should be able to tie these to surface 
pressure data to validate the noise source - somewhat 
similar in approach to that done in Ref. 6 tor trailing 
edge noise. We direct our attention to the edge pressure 
sensors on the suction and pressure sides. These would 
be the only sensors in the strong edge flow field and, at 
the same time, be in the near field of such a scattering 
phenomenon. They should therefore be representative 
of the source region. Note that the flap side-edge 
surface, between the suction and pressure sides, has 
generally lower velocity and its sensors (#1 through #9) 
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CFD Values (Simple Model Values) for M, 

0 = 0.17 



a =29° 



a =39° 


Sensor 

Number 

5 Inch 

M c 

p c < fc»g 

8 inch 

M c 

Pc deg- 

10 

.010 

.21 (. 222 ) 

77 ( 90 ) 

.020 

.22 (. 237 ) 

82 ( 90 ) 

11 

.010 

.21 .222 

90 ( 90 ) 
85 ( 90 ) 

.012 

.27 . 237 ) 

82 ( 90 ) 

12 

.010 

.23 (. 222 ) 

.025 

.29 (. 237 ) 

68 ( 90 ) 

13 

.030 

.22 (. 222 ) 

73 ( 90 ) 

.050 

.265 (. 237 ) 

80 ( 90 ) 

14 

.050 

.24 . 222 ) 

82 ( 90 ) 

.050 

.33 (. 237 ) 

76 ( 90 ) 

15 

.030 

.35 (. 285 ) 

54 ( 51 ) 

.045 

.32 (. 331 ) 

43 ( 46 ) 

16 

.040 

.33 (. 285 ) 

48 ( 51 ) 

.045 

.30 (. 331 ) 

46 ( 46 ) 

17 

.050 

.31 (. 285 ) 

52 ( 51 ) 

.060 

.28 . 331 ) 

60 ( 46 ) 

18 

.120 

.25 (. 285 ) 

i 58 ( 51 ) 

.150 

.24 (. 331 ) 

65 ( 46 ) 

30 

.010 

.13 (. 21 5 ) 

! 28 ( 34 ) 

.025 

.15 (. 253 ) 

10 24 ) 

31 

.020 

.19 (. 215 ) 

! 38 ( 34 ) 

.020 

.205 (. 253 ) 

28 ( 24 ) 

32 

.020 

.21 (. 215 ) 

35 ( 34 ) 

.020 

.205 (. 253 ) 

26 ( 24 ) 

33 

.025 

.22 (. 215 ) 

33 ( 34 ) 

.025 

.205 (. 253 ) 

25 ( 24 ) 

34 

.030 

.22 (. 215 ) 

32 34 ) 

.030 

.22 (. 253 ) 

25 ( 24 ) 

35 

.030 

.22 (. 215 ) 

31 ( 34 ) 

.030 

.21 (. 253 ) 

25 ( 24 ) 

36 

.030 

( 5 =. 0453 ) 

.20 (. 215 ) 

32 ( 34 ) 

.030 

( 8 - 0314 ) 

.205 (. 253 ) 

27 ( 24 ) 


TABLE 2. Calculated edge flow boundary-layer thickness and 
velocity values. 



Convective 

Velocity 


FIGURE 8. Flow above a surface sensor and an idealized 
shed instability wave geometry at edge of flap. 


are not considered here as representative of the source 
region (although they are in the near-field of such 
scatter). In Table 2, for the sensors indicated, values 
are given for the near boundary layer thickness 5, the 
corresponding Mach number M c , and flow angle /? ( , 
determined over planes parallel to the surface and at 
height z = S above the surface. The choice of 5 is 
partially subjective. It corresponds to the outer edge of 
the shear flow nearest the surface. An illustration of the 
velocity field U(z) above an edge sensor is shown in 
Fig. 8. The top view shown defines the angle /J c from 
the normal to the edge. The subscript designation c 
(for convective) is used to indicate the flow above the 
sensors is assumed to also represent any moving 
disturbance or flow structure that may cause noise- 
producing pressure scatter at the edge. The 
hypothesized convecting wake sheet illustrated in Fig. 8 
will be subsequently discussed. 

For both flap angle cases considered, Table 2 
indicates that <5, Af c , and /J ( remains generally 
invariant along much of the pressure side edge. The 
flow speed M c , as well as the cross-flow (or span wise) 


component of velocity M c Cosfi c , exceed the tunnel 
value of A/ 0 =.17 . On the aft (downstream of mid- 
chord) suction side edge, where the attached vortex 
flow comes off the surface past the edge, the flow 
velocities are even higher, reaching up to about twice 
the free-stream value. Forward on the suction side 
edge, the velocities are lower than those aft and the 
cross-flow diminishes greatly with flow skew angle p c 
approaching 90°. An unexpected result, to the present 
authors, for the CFD flow field is the lack of anticipated 
changes in 8 and M c values with changes in flap 
angle. Expected increases in M c did not occur with 
increased flap angle, even in regions further away from 
the surface. It should be mentioned that Ref. 28 noted 
that the solutions, while remarkably good overall in 
defining basic flow features, found disagreements with 
measured velocities on the order of 10 to 15%. 
Concerns about the thickness of the shear layer were 
also expressed. It was suggested in Ref. 28 that 
improvements may be needed with regard to grid 
resolution and turbulence modeling. Because of the 
importance of these parameters to the present effort, 
alternate calculations are made and are presented in the 
following section. The CFD solution, however, is 
utilized in providing a reference for primary flow-field 
features. 

Simplified edge flow calculations 

Simple aerodynamic modeling is used here to take 
into account Reynolds number and flap angle effects in 
the definition of boundary layer thickness and velocity 
values. This complements the description of the 
complex three-dimensional flow field given by the CFD 
solution. 

From thin airfoil theory 40 , the sectional lift per unit 
span L' equals 

c 

L' = pU 0 r = p£/ 0 J" ydx = q 0 cC N ( 1 ) 

0 

where p is the medium density and V is the airfoil 
circulation given an incoming stream velocity of (J 0 . 
The circulation density y of the vortex sheet defines 
the airfoil in the stream from the leading edge at x= 0 
to the trailing edge at c, where c is the chordlength (of 
the flap in the present case). The dynamic pressure is 
q 0 = p(J 0 / 2 and C N is the sectional lift coefficient 
defined by Eq. (1). Figure 9(a) shows a sketch of an 
inboard section of the flap where the flow is essentially 
two-dimensional. The velocity jump across the airfoil 
sheet is y(x) = u su - u pr , where u su is the velocity 
along the suction side and w is that along the pressure 
side. The mean or average velocity jump over the chord 
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(c) model flap vortex circulation model showing 
forward section cut A 1 and aft section cut B' 


FIGURE 9. Illustrations for simplified flap edge velocity 
calculations. 


is (u su -u pr ) mem = U 0 C N 12. The average values of 
u su and u pr are thus 

( w v u )mean = + C N f4) 

and (2) 

( u pr )mean = ~C N f 4) 


Along with the streamwise flow component above, 
the cross-flow component is required to estimate the 
edge flow. For this purpose. Figure 9(b) illustrates a 
somewhat heuristic flap vortex model, where the bound 
circulation of strength T of the flap is transitioned to a 
trailed wake vortex of the same strength at the edge. 
Figure 9(c) shows a forward section view (cut A’ in Fig. 
9(c)) where the vortex of core radius r 0 , is ’drawing* 
fluid at velocity u v from the pressure side into itself. It 
is assumed that the vortex does not greatly affect the 
edge flow on the upper suction side at this section. 
(Note that any contributions from ’’secondary” vortices 
are ignored.) For the aft section (cut B’) flow, the same 
basic vortex ’draws' fluid at the same velocity u v from 

both the pressure and suction sides. The maximum 
velocity of the vortex is at the core radius and equals 
m v . The velocity description chosen for the vortex is 
the Scully model 41 42 which distributes the circulation to 
simulate the effects of viscosity. The velocity due to 
the vortex is given in terms of the radius r from the 
center; u = (V 1 2/rr )[ r 2 /(r 2 + rj )]. (For an ideal point 


vortex, the bracketed term, [ ], would be unity.) 
Equating u v = u at r = r 0 , we obtain 

u v =U 0 C N cf*m 0 (3) 

where the strength F is obtained from Eq. (1). Using 
Eqs. (2) and (3), the velocities on the edge surfaces may 
be determined. For the pressure side, 


U pr=^j( u pr)Lan +U Y 

and (4) 

P = Tan~ ] l(u pr ) mea „/u v ] 


For the suction side, on the aft section where the vortex 
crosses to the upper surface, (J su is similarly expressed 
but the subscript su replaces pr. On the forward 
section, U su , however, is simply taken as (u su ) mean and 
/} =90°. In Table 2, values of these velocities (in terms 
of Mach number) and angles are given in parenthesis to 
compare with corresponding CFD values for the 
different sensors. The pr and su subscripts are 
dropped. The values were calculated using r 0 = 
0.4t max and previously mentioned value of C N = 1.213 
for a= 29° and C N = 1.567 for a= 39°. The maximum 
flap thickness is t max =0.55 inches. The value used for 
r 0 appears to give velocities and angles in nominal 
agreement with the CFD values, but unlike the CFD 
values have the physically expected flap angle 
dependence. 

Also listed in Table 2 (in parenthesis) are 
calculated approximate values of shear layer or 
boundary layer thicknesses 8 that one can compare to 
values that are determined from the CFD flow field 
previously discussed. The calculated values were 
determined by equations and extrapolated equations of 
Ref. 10, as defined below. The boundary layer 
thickness at the trailing edge for an untripped 
symmetric NACA 0012 airfoil at zero degrees angle of 
attack is empirically determined to be (Eq. (5) of Ref. 
10 ) 

_ c | q[ 1 .6569-0.9045/^/?, +0.0596{ Lt^R c ) 2 ] ^ j 

where R c =cU 0 /v is the Reynolds number based on 
chordlength c and v is the medium kinematic 
viscosity. At an angle of attack of a (taken as flap 
angle here) to the flow, the pressure side thickness 8 
can be related to 8 0 by 

5/8 q = io"° oi59 “ (6) 
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where a is angle in units of degrees. Note that Eq. (6) 
is newly determined here based on the data in Fig. 7 of 
Ref. 10. It replaces Eq. (8) of Ref. 10 in order to be 
more valid for large angles. The values of S using Eqs. 
(5) and (6), are 0.0453 in. and 0.0314 in. for 29° and 
39° at M 0 = . 1 7. These are somewhat larger than that 
listed in Table 2 for the sensors on the pressure and 
suction sides. 

In this paper, S is used as a normalizing parameter 
for surface pressure and overall flap noise spectra. Also 
used in normalizations are S values for roughened 
surfaces. This is approximated by the result of an 
interpolation of <5 0 between the untripped (Eq. (5)) and 
heavily tripped boundary layer cases of Ref. 10. The 
result is a replacement for Eq. (5) for lightly tripped 
surfaces, 

s = £ . ] o^ 1 787 — +0.0596( )* ] ^ 

Another normalizing parameter is an average M c 
value at the edge defined as 

‘j-'R +«? <8i 

a 0 


UNSTEADY SURFACE PRESSURES 
AND ANALYSIS 

Surface pressure spectra and acoustic source 
identification 

Figure 10 presents the unsteady surface pressure 
(auto-) spectra G s for four suction-side edge sensors 

and two pressure-side sensors (#32 and #34) for 
A/ 0 =.17. As previously stated, the data presented in 

this paper are limited to regions of flat frequency 
response for the sensors. It is a one-third-octave-band 
presentation, with the dB levels referenced to pi , 
where p 0 = 20pPa is the standard acoustic reference. 
In Fig. 10(a), for a= 29°, the levels are shown to be 
quite variable between sensors with the highest on the 
suction side at sensors #12 and #16. Referring to Figs. 
4 and 7, these two sensors are at opposite sides of the 
primary flap edge vortex on the suction surface. At 
a=39°, some relative level changes occur for all 
sensors. Referring to Table 2, one does not see any 
obvious correspondence between the velocity 
definitions over the surface at the sensors and the 
spectral levels. Figure 11 presents M 0 =.07, .11, and 

.17 data for two sensors normalized by q^.S I U c , where 
q c = pul 12. This type of normalization is common 


for surface pressure spectra under turbulent boundary 
layers, an example being Ref. 6. The values for U c are 
determined as f/^and U s from Eq. (4), as was done to 

obtain values for Table 2. The values for S are 
obtained from Eqs. (5) and (6). These values depend on 
velocity and flap angle, but not chordwise location. It 
is seen (by the degree of data coalescence) that while 
tunnel velocity dependence is partially captured, it is 
not consistent between flap settings. The spectral 
shapes for each sensor apparently depend greatly on the 
particular flow phenomena occurring about it - 
therefore, the local flow phenomena apparently change 
with angle and velocity variations. The parameter 
grouping does not capture this. It is noted that the use of 
the S and U c values from Table 2 based on the CFD 
results produces no improvement in normalization 
success. It is expected that any such normalizations 
should be more successful at sensor positions more 
inboard, away from the edge. 

(a) a =29° 




FIGURE 10. Surface pressure spectra for pressure sensors at 
edge of flap for M 0 - .1 7 for different flap angles. 

Figure 12 shows some spectral characteristics that 
can provide a basis for a noise mechanism hypothesis. 
The auto-spectra G s are given for edge and inboard 
sensors for both the suction and pressure sides. The 
spectral resolution is Af= 244 Hz but the levels are 
referenced to a A/ = 1 Hz bandwidth. For the suction 
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1 io r 






FIGURE 1 1 . Normalized pressure spectra for different angles 
and flow speeds. 


Sensor 




FIGURE 12. Pressure spectra and cross-spectral phase 
relationships for edge and inboard sensors on suction and 
pressure sides of flap. The spectral resolution is A f - 174 Hz, 
but levels are referenced to A/ = 1 Hz. 

side, the inboard sensors #21 and #23 are comparable in 
level to one another and are lower than the edge sensor 
#16 levels by about 15 dB at 5 kHz. For the pressure 
side, the levels of the inboard sensors #40 and #42 are 
comparable, to one another, and are lower than edge 
sensor #34 by over 15 dB at 5 kHz. The levels for the 
farther inboard sensor #46 are even lower. This 
characteristic of increased surface pressure spectral 
levels, as the edge is approached from inboard is 
counter to that found for the classical turbulent- 
boundary-layer (TBL) trailing-edge (TE) noise scatter 
problem 6 . (There, of course, the radiating edge is the 
trailing edge rather than the present flap side edge. In 
that case, the hydrodynamic (TBL) pressure field is 
intense upstream of the edge. But very near the edge, 
the levels decrease, because of pressure scatter (near- 
field noise) that prevents a pressure differential at the 
edge.) Therefore, the noise level behavior of Fig. 12 
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suggests a different mechanism than that for the TBL- 
TE noise problem. This is correspondingly true for the 
phase behavior to be discussed below. 


u h 



FIGURE 13. Hypothesized flow and acoustic features affecting 
unsteady surface pressures at sensors. 

The hypothesized mechanism for the present flap 
edge problem is illustrated in Fig. 13, which shows 
flow-field influences on unsteady surface pressures. 
The subject sensors of Fig. 12 are shown mounted in a 
"section cut", the edge sensors are located at y A and the 
others are at y B and y c of Fig. 4. The edge #16 and 
#34 sensors are not aligned chordwise. The 
conceptional illustration is consistent with the shear- 
layer instability models^ ’ for noise production. 
Shear layer instabilities are shown being shed at 
velocity U c at the edges near sensors #16 and #34. The 
velocities U c , U hJ and U' h should be of the same order 
of magnitude to that of velocity u v of Fig. 9(c). The 
flow is relatively smooth on the pressure side but 
becomes turbulent, and convects at say velocity U h , as 
it moves around the edge towards the suction side. The 
suction side turbulence is within the fringes of the 
primary vortex and convects at velocity U' h above the 
surface. Figure 13 represents the shear layer instability 
noise source as dipoles (that are distributed chordwise). 
A portion of the noise that is radiated travels along the 
surface in the spanwise direction at the speed of sound 
a 0 . Both dipoles radiate to both sides with opposite 
signs (180° out of phase). Of course Fig. 13 is a 
sectional presentation. In reality there are a number of 
independent dipoles radiating at different sections - the 
effective number of which depend on disturbance 
correlation scales. 

The observation that the noise levels of Fig. 12 are 
diminished as distance from the edge is increased is 
consistent with the model of Fig. 13. The phase 
behavior is also consistent as is now shown. The phase 
(p between sensor #34 and the inboard sensors, on the 
same pressure side, is normalized by subtracting 
(fd/a 0 )‘ 360, where / is the frequency and d is the 


distance from #34 to the other sensors (that one 
determines from Table 1). A result of zero degrees 
would show that the correlated components of the 
respective sensors are the same signals that are simply 
time-delayed at the speed of sound. The results shown 
in Fig. 12 show this to be generally true, except for a 
30° to 60° offset. However, note that in Ref. 6, the 
scatter term for TE noise was found to have a 45° offset 
in scatter-pressure phase (in addition to that due to time 
delay) between the edge near-field and a point away 
from the edge. Hence, one would expect a similar 
functional form for this scatter problem to that of Ref. 
6. Additional confirmation of the conceptional model 
is provided in Fig. 12 by the phase of sensors on 
opposite sides of the flap. The phases are not 
normalized. It is seen that for sensors away from the 
edge, there is roughly a 180° shift over much of 
frequency range. At the edge sensors #16 and #34, an 
approximately 180° phase is attained near 4 kHz. At 
lower frequencies, the phase is dominated by 
hydrodynamic convective effects, symbolized in Fig. 13 
by turbulence moving at U h . The linear phase slope, 
starting near zero frequency, suggest that U h between 
the sensors is about .22 times a 0 . Similar values were 
found for U' h using the sensors on the suction side 
surface. This compares favorably with computed (CFD 
and simplified) convective velocities. It was found that 
hydrodynamic-convection-effects generally dominate 
the phase relations between most edge sensors, as well 
as those over the suction surface. Most of these effects 
are not directly related to noise production. However, 
this does not mean that the individual sensor auto- 
spectrum is not dominated by noise related effects - it 
just means that the hydrodynamic effects are larger 
scale and thus correlate better over distance. The 
present edge sensors are deemed too far apart for cross- 
spectra to determine pertinent noise source information, 
such as scale lengths. 

Figure 14 serves to summarize key features of the 
surface pressures. The figure presents chordwise 
distributions of integrated surface pressure levels for 
sensors, at sections A, B, and C of Fig. 4, for both 
pressure and suction sides of the flat edge flap for 
a =29° and 39°. The levels result from integrating the 
spectra from 4.0 to 13.5 kHz. The low frequency limit 
of 4.0 kHz was chosen in order to de-emphasize purely 
hydrodynamic effects (see discussion above). For both 
flap angles, the levels at the suction side edge (A) 
appear peaked in the general locations of the chordwise 
extremities of the vortex on the suction side. On the 
pressure side, one peak is observed near the flap mid- 
chord. The inboard sections (B) and (C) have levels 
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FLAT EDGE 



depicted in Fig. 13. The higher levels on the suction 
side for the inboard sensors, with respect to the edge 
levels, are expected to be due to the presence of strong 
turbulence on the suction side. 

Figure 15 presents chordwise level distributions for 
the round edge flap sketched in Fig. 2. The character of 
the distributions is somewhat similar to that of the flat 
edge, but the section (A) sensor levels are reduced to 
nearly that of the inboard sensor levels. This is because 
they are further inboard of the edge than was case for 
those on the flat edge flap. The noise mechanism 
details that are depicted in Fig. 13 are not directly 
applicable to this round edge flap case. However. How 
shear / boundary layer instabilities still should be the 
basic mechanism, but with a different scattering 
geometry and correlated fluctuation length scales. 

Coherent Output Power analysis 


FIGURE 14. Chordwise distribution over FLAT edge flap 
region of band-limited overall surface pressure levels for 
a =29° and a =39°. A/ 0 = .17. 


ROUND EDGE 



x/c 


x/c 


FIGURE 15. Chordwise distribution over ROUND edge flap 
region of band-limited overall surface pressure levels for 
a =29° and a =39°. A/ 0 =.17. 


that are relatively uniform over the chord and 
substantially lower than the edge levels. Section (C) 
levels are generally lower than those of section (B). 
This is consistent with the noise mechanism modeling 


The surface pressure levels in Figs. 14 and 15 do 
not necessarily indicate the local noise-source strength 
distribution. The inboard levels contain substantial 
contributions from hydrodynamic convection effects 
and noise radiating along the surface * not sources of 
noise. (Note that strictly speaking, in terms of the 
Ffowcs Williams and Hawkings equation 43 , the inboard 
pressures are by definition part of the noise region that 
should be accounted for. This is discussed in the next 
section.) The edge sensors for the Hat edge flap, 
however, are in the near-field of the source and can 
represent the source, except to the extent that it is 
affected by non-radiating fluctuating pressure 
components. This section is concerned with providing 
some measure of the noise source distribution along the 
flap edge. 

The Coherent Output Power spectrum 44 is defined 
as 

CO pJ^- = r l s G s (9) 

where, for present purposes, this is the spectrum of the 
surface pressure sensor (subscript s ) output that is 
coherent with the SADA array (subscript a) output 
when steered to this surface sensor(s) on the flap. The 
auto-spectra are G a and G s . The cross-spectrum 

between signals is G a s and the coherence is y~ s . The 
phase associated with COP s is the cross-spectral phase 
(p a s . In the present data processing, a data-block time- 

shifting procedure is used to avoid serious bias and 
statistical errors, as well as to put the data in a useful 
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phase format. The microphone raw time data are time 
shifted (offset) by an amount close to the value of T as , 

which is the time required for noise to be transmitted 
from the sensor to the array. The array processing 
shear-layer refraction correction code determines t a v . 

(The T a s values were evaluated to be accurate within 

less than the time it takes to acoustically travel .25 
inches.) Final adjustment to obtain the full T a s shift 

effect, for each individual microphone and sensor, is 
done in the frequency domain. With regard to the cross- 
spectrum, this effectively puts the source region in 
’’retarded" coordinates, where the phase related to noise 
transmission time is removed, that is 

(C U .A„., =G a y jm °-< 

and ( 1 0) 

(<Pa.s) r„, =9a.s-<Wa.s 

where the radian frequency (0 = 27tf . 


FLAT EDGE 




FIGURE 16. Coherent Output Power spectral processing 
results, for the flat edge flap, relating the sensor #34 pressure 
measurement and the SADA noise measurement for A/ 0 = .17 
and a- 39°. The spectral resolution is A f - 34.88 Hz, but 
levels are referenced to Af= 1 Hz. 


FIGURE 17. Chordwise distribution over FLAT edge flap of 
band-limited overall surface pressure levels and related COP 
levels for or =29° and a=39°. Af 0 = .17. 


ROUND EDGE 



x/c x/c 


FIGURE 18. Chordwise distribution over ROUND edge flap of 
band-limited overall surface pressure levels and related COP 
levels for or =29° and a =39°. A/ 0 =.17. 
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Figure 16 shows results of COP processing for 
sensor #34 (pressure side) with respect to the output of 
the SADA when steered to the sensor for the test case 
listed. The SADA is positioned at 0 = 107°. The two 
auto-spectra, the cross-spectrum, and the COP s spectra 
are shown. (Note that the differences in surface 
pressure auto-spectrum smoothness are related only to 
an application of a calibration transfer-function. As 
previously mentioned, the surface pressure data should 
be accurate below about 13.5 kHz.) The difference in 
levels between G s and COP s is 10 Log{y^ x ), which is 
shown in the figure. As is shown in the next section, if 
a surface element (represented here by sensor #34) is a 
direct radiator to the noise field, one would expect the 
phase (<p a s ) T to equal a constant -90°. The phase is 

shown in Fig. 16 to be generally constant, but at a phase 
value on the order of -40° or -50°. However, this is 
consistent with a 45° offset expected for an edge in a 
scatter field as mentioned in the last section. This 
phase behavior will be again discussed. 

Figure 17 presents chordwise distributions of 
integrated surface-pressure COP levels for the edge 
sensors (along section A). For comparison, the auto- 
spectra distribution from Fig. 14, for the same 
integration frequency range, is shown. No data from 
inboard sensors are shown because these are assumed 
not to be in the noise source region. The COP 
distributions represent more realistic distributions of 
noise source strength distribution than does the auto- 
spectra. The COP results eliminate that portion of 
surface pressure that is not related to the noise field. So 
non-radiating hydrodynamic fluctuation contributions 
are removed, which lowers the COP levels with respect 
to the auto-spectra levels. Also, however, there is an 
additional cause of the lower levels for the COP results. 
The correlation area that each sensor represents is small 
and there are a large number of correlation areas over 
the flap edge which contribute to the total noise. For 
example, if there were uniform noise source strength 
(say with no non-radiating hydrodynamic fluctuations 
present) and the COP were uniformly 20 dB lower than 
the auto-spectra, one could hypothesize that there are 
"effectively'* 100 independently radiating noise source 
areas across the flap edge. Figure 17 shows two COP 
distributions. One is where phase is not considered in 
the integration (straight pressure-squared type 
summing) of the frequency bands. The other 
distribution is where the phase is used in the integration 
(vectorial type summing). The latter is lower in level 
and is the preferred presentation because those portions 
of hydrodynamic and/or acoustic fluctuations, which 
are related to noise production only in an indirect way. 


are substantially eliminated. Consider the COP peak 
near 20% chord for the phase-suppressed distribution 
on the suction side for a =39°. Phase data (not shown), 
indicate that the fluctuations which cause the high COP 
levels were related to turbulence and/or noise, that are 
in turn correlated with noise production at another 
portion of the edge. Since the phase «p a v ) r ^ has to be 
constant for that portion of COP related to the direct 
radiation from the sensor to the microphones, the 
summing of COP bands vectorially can substantially 
remove (bias against) the "indirect" contributions to 
COP. This should make the COP more representative 
of the actual source distribution. The COP distributions 
in Fig. 17 show that, for both of the flap angles, the 
noise is most strongly radiated near 65% chord on the 
suction side and near 50% on the pressure side of the 
flap edge. 

Figure 18 shows COP results for the round flap 
edge, in the format of Fig. 17. The same comments 
apply here, as for the flat edge flap, except that the edge 
sensors may be not be fully in the source region, as 
mentioned for the auto-spectra integrated level plots of 
Fig. 15. Still, one can make the statement that the 
chordwise noise source appears clearly and strongly 
located near 60% chord for both pressure and suction 
sides for both flap angles. 

NOISE PREDICTIONS BASED ON SURFACE 
PRESSURE SPECTRA 


Causality spectra prediction 

A causality spectral approach is developed in this 
section that helps establish the relative and quantitative 
importance of the noise source regions of the flap. 
Previous success has been found by Siddon 4 \ using 
cross-correlation methods and the acoustic relationships 
of Curie 46 , in determining surface noise source 
distributions. This causality approach employed cross- 
correlations between surface pressure sensors and 
microphones. For several simple and small surface 
shape cases under Siddon’s study, the method provided 
a physical characterization of the noise source. 
However, when the sources were non-compact, 
acoustically or aerodynamically, the phase variations 
greatly hindered correlation function interpretation and 
their usefulness. In this paper, we revisit the causality 
idea using spectral methods and for the first time 
validate a causality prediction with measurement tor a 
distributed source. 
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The noise field is given by the Ffowcs Williams 
and Hawkings equation 43 (a form similar to Curie’s 
equation but generalized for arbitrary fluid and surface 
motion). For low-Mach number flows and for surfaces 
with steady (or no) motion with respect to the observer, 
distributed volume quadrupole and surface monopole 
source components in the equation are negligible. 
Assuming surface shear stresses are small compared to 
local surface pressures, the following equation form can 
be found, relating the acoustic pressure p a (x,t) at 
location x and time rto the surface pressure p s (y, T) at 
position y and retarded time T 


A, (*.*) = 


1 

4 7ta 0 r 


r nr fr 

dp A y.T)' 

J (l-M 0 r/r) 2 

dr 


S(y) 


( 10 ) 


Figure 19 shows the geometry of the flap in the 
open jet tunnel flow of Mach number M 0 = U 0 /a 0 , 
where U 0 is the tunnel test section velocity. The 
elemental surface area dS( y), at y with normal n, is 
seen with respect to a ray path of length r = |r|, where 
r = x-y. Shear layer corrections performed in the 


Pa(**t) = 


4 m 0 rD 2 


i 


CosQ\ 


S( y) 


[ ^v(y.T) ' 

dx 


L 


\dS(y) 


(13) 


where D=I-M 0 -r/r = l-M 0 x/|x| and 

CosO — n • r / r . The retarded time becomes 


r' 

r = (/ ) (14) 

«o 

where r' ~ rD 2 , which is shown in Fig. 19. The 
distance r', which is used in the shear layer processing 
code 35 of this study, is that between the effective 
observer position and source position in an ideal 
quiescent acoustic field. In this coordinate system, 
CosO is replaced by CosO ' = n * r7 r f . Upon defining 
the Fourier transform 

3( P s (y, T)] = J p s (y, T )e~ jB *dt = PJy,co) (15) 

PAV' 1 ~ r ' t a (>)i = -j(oP s (y,w)e~ ja>r ' !an (16) 


Mean 



so the Fourier transform of the acoustic field is 


P a (x,(o) - — f CosQ’P s { y, co)e- Jkr dS( y) ( 1 7) 

4^a n r J 

° Si y) 

where k = 0)f a 0 . Multiplying both sides of Eq. (17) by 
the complex conjugate of P a (x,(o) and taking the 
ensemble average, we get the auto-spectrum of the 
noise field 


FIGURE 19. Flap surface geometry with respect to ray path to 
noise observer. 

array processing code corrects the results at the actual 
SADA out-of-flow position to an ’effective' observer 
position x in an extended flow field without a shear 
layer. The medium speed of sound is denoted by a 0 . 
The retarded time is defined implicitly by 

r = ;-|r-U 0 r|/a 0 (12) 

as pointed out by Guo 33 . In the far-field limit for x, 


G a ={P ll *(x,co)P a (x,G))) 

= j~p J Cos6\P a \x,(D)P s ( y, (0))e~ jkr 'dS(y) ( 1 8) 

° r Sty) 

Now upon evaluating this with respect to elementary 
surface areas AS(y), the following results 

C„=(«) 

(19) 
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The components, making up the sum in the right hand 
side of Eq. (19), are defined here as the causality 
spectra associated with the surface area elements AS(y) 
and represent their contribution to the total noise G a . 
The auto-spectrum G a must be a positive real quantity, 
while the components are complex quantities. Equation 
(19) is valid under the present assumptions, in 
particular that only surface dipole noise sources are 
important, and that the sum is taken over infinitesimal 
AS(y) areas over all of S(y). Equation (19) includes 
all the hydrodynamic pressure fluctuations (correlated 
or not between different A5(y) areas), as well as related 
pressure scatter (near field noise). It is important to 
note that non-radiating turbulence hydrodynamic 
effects, with their various complex phase contributions, 
may greatly dominate the individual terms. The full 
evaluation of the sums are needed to completely self- 
cancel the non-radiating contributions. Such concerns 
not only apply to Eq. (19) but also to the Ffowcs 
Williams and Hawkings equation in general. 

In those practical application cases where AS(y) 
must be finite in size and limited in number, the validity 
of Eq. (19) requires that AS(y) be chosen and 
interpreted carefully. As previously noted in 
discussions pertaining to Fig. 17, many of the pressure 
sensors of the present problem were found to be indeed 
dominated by non-noise producing hydrodynamic and 
acoustic effects, for which portions are correlated with 
the noise. The use of the edge sensors only for the flat 
edge flap should reduce extraneous 'noise' in any 
prediction and data comparison. The model then 
reduces to a line of dipoles along each side edge, with 
the presence of the inboard flap surfaces not included in 
the solution for the acoustic radiation. To account for 
the "half-baffle" acoustic effect of the inboard surfaces, 
the solution can be multiplied by a factor of two (2). 
This should give an accurate presentation for small 0 ' . 

For the present study, Eq. (19) is evaluated using 
data from the flat-edge flap. Our attention is restricted 
to the edge sensors #1-18 and #30-36 as representing 
the source region. The following relation is evaluated 

C„ -£<«), 

i 

. 7T 

= - jgLy [ c °sO’* _L ]/(py kr ' )• p\u I (20) 

4 7ta {) j r ' 1 

where the additional factor of 2 included is discussed 
above. The cross-spectral term containing the retarded 


time shift is identified as (G u v ) r ^ of Eq. (10). The 
area AS(y) has been replaced by Lf } . We take L to be 
the sum of correlation length scales that sensor i 
represent and, nominally, the sum of all L equals the 
edge circumference. The correlation length scale in 

the flap longitudinal direction (spanwise direction in 
Fig. 4) is taken as U c /rjco. This relationship for scale 

length is often used in turbulent boundary layer (TBL) 
pressure scaling. In the present study, the available 
data provided no satisfactory means to determine values 
of 77. In Ref. 6, r\ values for a TBL was found to be 
.14 and .19 for the cases examined. For the following 
predictions, we use the value rj =.3. This choice is 
discussed in a following section. 

Causality predictions and comparisons with 
measured noise is presented in Figure 20 for the flat- 
edge flap at M 0 = .07, .11, and .17 with cx — 29°. The 
SADA is positioned at 107°. Two causality prediction 
results are shown. The predictions based on Eq. (20) 
are the prediction curves showing the lower levels. The 
associated phase for the A / 0 =. 1 7 speed case is shown 
plotted below the spectra. Ideally, the phase of G a 
should be zero. However, the phase is seen to be 
approximately 30° to 50°. This is consistent with the 
aforementioned 45 s phase shift in the edge sensor 
pressures located in the near-Field of the edge scatter. If 
G a were fully evaluated over the whole surface (not 
possible with limited data) then the phase would be 
expected to be near zero. The other causality prediction 
results shown are where the phases in Eq. (20) are 
suppressed for each individual contribution. Forcing 
phase to be zero removes additive random phase error 
in the cross-spectra and errors related to time delay 
variability for each sensor. It may, however, add bias 
error of an unknown amount by adding correlated 
components, when they would otherwise properly 
cancel. Still the results appear to compare well over the 
whole spectra range. For the a = 39° case, the same 
prediction comparisons are made. This is shown in Fig. 
21 in the same format as Fig. 20. For this case, the 
predicted levels are lower than measured levels, 
particularly for the causality predictions where phase is 
included. Note correspondingly that substantially 
larger phase variations are seen compared to those in 
Fig. 20. This indicates that these predictions based on 
only the edge sensors are less representative of the total 
noise production. The influence of the burst vortex, 
associated with this a= 39° flap angle case, may cause 
the noise source region to be more distributed over the 
surface. 
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Frequency (kHz) 


FIGURE 20 Causality noise prediction comparisons for flat 
edge at different tunnel speeds for a =29°. T) =.3, Af =244 Hz. 



Frequency (kHz) 


FIGURE 21. Causality noise prediction comparisons for flat 
edge at different tunnel speeds for a =39°. 77 =.3, Af =244 Hz. 
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Scatter edge noise prediction 

Brooks and Hodgson 6 used the trailing edge noise 
theory of Howe 47 and measured surface pressures to 
predict trailing edge noise due to the passage at the 
edge of turbulent boundary layer flow. The present 
edge noise problem is that of shear layer wake 
instability and resultant shedding of unsteady vorticity 
from the edges. Howe's theory was developed 
primarily with the former case in mind. However, the 
solutions should be generally valid here as long as one 
restricts attention to pressures in the immediate vicinity 
of the edge and that certain parameters can be properly 
defined. Figure 8 is used to illustrate several parameter 
choices in doing this. For the velocity field U(z) 
above an edge sensor, a local maximum of magnitude 
U c is reached at a height <5. These values correspond 
to those listed in Table 2, along with corresponding 
skew angle fl c . One hydrodynamic wavenumber 
component of the instability wake sheet is shown in 
Fig. 8 as being shed from the edge region. It is 
assumed that this sheet perturbation (shown as a 
corrugation) convects at this same speed U c and angle 
fi c after it leaves the surface. It is also assumed that the 
edge sensor is in the very near field of the edge 
shedding and that the local edge thickness is much 
smaller than the related acoustic wavelength. 

Equation (72) from Howe 47 (Eq. (32) from Ref. 6) 
gives the noise spectrum G a at a location in the far 
field due to the trailing edge (TE) noise from a thin 
plate of length L, in terms of the TE pressure field. 
This is, in the present terminology, 


2 ( |Co.Tg[5w 2 (;»/2)Cfty/j, 

a 0 J(1 + M„ r ) 2 (l - Mv Jd-M^lCosdl) 

J n , c aSin 6/ a 0 , (o)d}J x (21) 


where 0 is the observer angle defined in Fig. (19). The 
observer azimuth angle & is measured from the 
negative of the y axis shown in Fig. 4 ( & equals 180° 
along the spanwise surface and equals 90° normal to 
and above the flap surface). The Mach number terms 
are M 0r = U Q Sin$/a 0 , M Ur = U c Cos& I a 0 , and 

M ^ = U c Cos/3 c / a 0 . Equation (21) ignores one of the 

Mach number terms of Ref. 47, see Ref. 6. The integral 
is taken over a pressure wavenumber spectrum function 
U te with respect to the wavenumber for the 
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direction normal to the edge. In evaluating the integral, 
we use the definition of fl f( , from Ref. 6 (Ret. 47 uses 

an alternate but consistent definition), to obtain 

Jr (22) 

— oo 


where (G s ) te is the surface pressure spectrum at the 
edge and is the correlation length scale in the lateral 
(edgewise along chord) direction. The assumed form of 
is 




U c Cosf} ( . 

geo 


(23) 



(a) a =29° 


The noise is predicted as a sum of contributions 
from individual edge lengths each represented by an 
edge sensor. Similar to Eq. 20 for the causality 
prediction, the total noise is 

(24) 

i 

where each [GJ, is determined using Eqs. (21), (22), 
and (23). As mentioned for the longitudinal scale factor 
in the last section, correlation scales for the present 
mechanism were not determinable from the present 
data. In the following predictions the value of g = 2.0 
is chosen compared to a value of 0.6 (measured for the 
problem of a TBL pressure field) used in Ref. 6. This 
choice and other assumptions are discussed in the next 
section. 

The predictions are compared to measured noise 
for the two Hap angles at different tunnel speeds in Fig. 
22. The measured noise spectra are the same as 
presented in Figs. 20 and 21. The chosen value of g 
results in good agreement for the 29° flap angle, as well 
as the low speed conditions for the flap angle of 39°. 
The higher speeds for 39° show predictions to be lower 
than measured. This general trend, of course, was also 
found for the causality prediction comparisons using 
the same sensors. One can then suggest that for the 39° 
flap case, where the vortex is known to burst, the edge 
sensors may not fully represent the noise production 
region. 

Discussion of prediction results 

The present predictions strongly support the basic 
noise mechanism model of Fig. 13 for the flat-edge flap. 



(b) a =39° 


FIGURE 22. Scatter-theory noise prediction comparisons for 
flat edge flap at different tunnel speeds for two flap angles. 
g = 2.0. A/ =244 Hz. 

The good prediction results show that the edge sensors 
can properly represent the noise source region. 
However, the prediction and data comparisons suggest 
this may not be fully true for the vortex bursting case, 
where the noise is under-predicted. It is likely that the 
source is just more distributed over the surface. (At the 
low speeds of the present test, it is unlikely that any 
volume noise source terms contribute significantly 
even with vortex breakdown.) The two predictions are 
different in their input requirements. The causality 
prediction and the COP analysis can be regarded as 
being related, because of their dependence on cross- 
spectral processing between the noise and the surface 
pressure sensors. Note that the causality predictions 
require little flow information input, with velocity 
entering only through the definition of correlation 
length. The causality prediction is primarily one of a 
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noise-correlated force distribution definition, modeled 
here as a line-dipole taking the surface into account 
acoustically only through the multiplicative factor of 
two. The scatter edge noise prediction is more of a 
"full" prediction, which more properly includes source 
directivity in the solution. It requires more flow 
information than the causality prediction. However, for 
both predictions, the lack of knowledge about pertinent 
correlation scale lengths required assumptions. The 
values for decay factors Tj and g were chosen to render 
good overall quantitative comparisons. These 
correspond to scale lengths of t x = .4 in. and *=.06 
in. at 5kHz. The ratio ^ = 7 is compared to a value 

of about 4 for the different scatter noise problem (TBL- 
TE) of Ref. 6. Still, the ratio IS *1.5 appears to be 
’’reasonable" for the present mechanism, although 
/<5 = 10 may not be. Subsequent investigations 
should reexamine the correlation scale issue, 
particularly key parameters such as the disturbance 
velocity magnitude U c and skew angle Still, 

uncertainty about correlation length presentations does 
not undermine the basic physical understanding and 
theoretical context gained from the present prediction 
comparisons. 

NOISE DIRECTIVITY AND 
SPECTRA SCALING 

Noise source distribution from array scanning 

Acoustic results from the array are shown in Fig. 

23 for the flat and round flap edges. The results are 
obtained from the SADA by electronically scanning a 
plane projected through the airfoil main element 
chordline. The position of the SADA corresponds to 
the model being in an ’’over-flight" position. An outline 
of the main element is shown with the leading edge at 

24 in. and trailing edge at almost 40 in. in tunnel 
coordinates. The flap is on the right and the edge is 
seen centered in the picture. The dB levels shown are 
the outputs of SADA when it is steered to the scanning 
locations. The contour levels are highest at the flap 
edge location. These levels at the flap edge have been 
shown 36 to be the levels that a single microphone would 
measure from the flap edge. The rapid roll-off in levels 
away from the flap edge shows the sharpness of the 
array in rejecting unwanted extraneous noise from 
regions other than the edge. The contours shown are 
for 12.5 kHz one-third octave levels. Because of the 
microphone shading algorithm methodology 35,36 , other 
frequencies from 10 to 40 kHz show similar spatially- 
invariant patterns. At lower frequencies, the resolution 
decreases (patterns widen) and the array rejection of 
extraneous noise is reduced. Still the levels from the 


vicinity of the flap edge should have little 
contamination for frequencies above approximately 3 
kHz. The spectral output of the SADA should hence 
represent only that noise which is radiated from the 
flap-edge region. Noise directivity (shown in the 
following section) is mapped by placing the SADA at a 
series of elevation and azimuthal angles, while 
maintaining a constant distance of 5 ft. from the flap 
edge region. 


f 1/3 = 12.5 kHz 



FIGURE 23. Noise source distribution contours over the flap- 
edge region using the SADA for flat and round edges at the two 
flap angles. SADA position is 0=107° and 0=0°. One-third 
octave levels for f m = 12.5 kHz. 

Directivity 

Figure 24 shows the model with the flap-edge 
directivity contour mapped over a spherical surface, 
defined by the SADA positions. The measurements are 
for the flat edge flap model for a =39° and A/ 0 =0.17. 
For the 6.3 kHz one-third octave frequency band 
shown, the directivity on the pressure side of the model 
is most intense “underneath' 1 the model. This is the 
side that an observer would "see” when an aircraft flies 
overhead. On the suction side of the model, the levels 
are less but are seen to increase in the downstream 
direction. Figure 25 are pressure-side directivity maps 
for a= 29° and 39° and selected frequencies ranging 
from 3.2kHz to 40 kHz. These maps are flattened 
versions of the spherical surfaces shown in Fig. 24. 
The positive azimuthal angles \| / are on the flap side of 
the model. The elevation angles <J) with the smaller 
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values (at the top of the plots) are in the 
downstream direction. For this flat flap-edge 
configuration, the directivities have a simple dipole-like 
shape at lower frequencies with the dipole axis oriented 
inline with angle <p at 90° to 107° and \|/=0°. For 
higher frequencies (>12.5 kHz), the directivity has a 
more “baffled” dipole character with stronger levels on 
the flap side (\j / negative). The directivity for the flange 
edge is shown in Fig. 26, where at lower frequencies 
the levels and patterns are similar to the flat edge 
results. However, at higher frequencies, the directivity 
patterns are somewhat more complicated and the levels 
are higher. Results for the round edge flap are even 
more complicated with generally higher levels and 
multiple directional peaks; suggesting a more complex 
edge noise source (pressure scatter region of surface) 
than that found for the flat edge. With the application 
of grit to the pressure side of the round edge flap, the 
levels decrease to approximately those seen for the flat 
edge flap. However as Fig. 28 shows, the directivity 
remains just as complex as the round edge without grit. 



FIGURE 24. Flap edge noise directivity over 3D "surface" 
defined by the SADA measurements. One-third octave levels 
for / 1/3 = 6.3 kHz. 


Spectra and scaling 

One-third octave spectra for the three flap 
configurations are shown in Fig. 29 for the SADA 
located at 0=107° and 0=0°. The round edge is seen 
to be the loudest configuration at low frequencies but is 
the quieter at higher frequencies than both the flat and 
flanged edge flap. For the flanged edge flap, the broad 
spectral peak at higher frequencies is particularly strong 
(and thus troublesome). Figure 30 shows spectra for 
different tunnel speeds for all three configurations, plus 
the flat and round edge configurations with grit applied. 
The spectral level dependence on tunnel speed is seen 
to be the most important, followed by flap angle and 
then by the application of grit. The grit serves to trip the 
boundary layers causing them (and thus the off-the- 
edge shear layers) to become more turbulent and 
thicker. Noise levels are reduced. 

The spectra of Fig. 30 are scaled by normalizing 
the levels and frequencies with M Cavc from Eq. (8) and 

8 from Eqs. (5) and (6) for the basic configurations, 
but from Eqs. (7) and (6) where grit is applied. Figure 
31 presents this scaling. The levels are referenced to 
the fifth power of the Mach number term. The levels 
are not taken to depend on <5. This normalization thus 
ignores a slight decrease in low frequency noise found 
for the flat edge when grit is used, but is consistent with 
negligible change in low frequency noise for the round 
edge when grit is used (see Fig. 30). The primary effect 
of thickness 8 is (taken as) to simply shift the noise 
levels to a lower frequency based on the Strouhal 
number f }/ *8/U r . Note that the normalization 

brings the spectra for the flat edge, with and without 
grit, into good general agreement. However, there 
appears to be some speed or Reynolds number 
dependence in spectral shape not accounted for. The 
normalization for the round edge is quite successful. 
The spectra data for the with- and without-grit cases 
appear well matched and coalesced. The spectral 
normalization for the flanged edge is also generally 
good. A lack of coalescence is seen over the broad 
high-frequency peak, likely related to the flange cavity. 
For all configurations, significant success is found in 
capturing the flap angle a dependence through the use 
of velocity U ( - AVC ■ Because this velocity depends on the 

flap C N , it provides the appealing connection between 
noise and flap loading. 
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V (deg) y (deg) 


FIGURE 25. Flap edge noise directivity over projected surface 
for FLAT edge flap. M 0 = .17 and two flap angles ( a =29° and 
39°) for different one-third octave frequencies. 



V (deg) y (deg) 


FIGURE 26. Flap edge noise directivity for FLANGE edge flap 
for conditions of Fig. 25. 
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ROUND EDGE 


ROUND EDGE W/ GRIT 
29° 39° 29° 39° 




FIGURE 27. Flap edge noise directivity for ROUND edge flap FIGURE 28. Flap edge noise directivity for ROUND edge flap 
for conditions of Fig. 25. with GRIT on surface for conditions of Fig. 25. 
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FIGURE 29. Noise spectra for the FLAT, FLANGE, and 
ROUND edges for two flap angles at Af 0 = 17. 



CONCLUSIONS 

The flat (or blunt), flanged, and round geometries 
are flap edge configurations in actual use on aircraft 
today. Spectra and directivity are presented for each to 
form the basis for semi-empirical predictions. New 
scaling methods are given. Significantly, with some 
exceptions noted, the noise spectra for each edge 
configuration are very successfully scaled for flow 
speed, flap angles, and surface roughness, using the 
simplified flow and boundary layer calculation methods 
developed herein. 

This study also provides an experimental and 
theoretical validation that, for the flat flap edge, shear 
layer instability and related pressure scatter is the 
dominant noise mechanism. For a higher flap angle, 
measured noise levels exceed predictions, which 
suggest additional contributions from surface sources 
that are not localized to just the immediate edge region. 
Such a source is likely related to the vortex bursting 
that occurs at the higher angle. As with the flat edge, 
shear-layer instability effects should be responsible for 
noise for the rounded flap edge, but the geometric 
features for the mechanism are different, producing 
different pressure scatter patterns and noise directivity. 
The flanged edge flap has an additional (and 
troublesome) high frequency noise contribution, likely 
due to cavity-type effects. 

Several newly applied diagnostic tools based on 
cross-spectra are used to determine the character of the 


(a) FLAT 



FIGURE 30. Flap edge noise spectra for (a) FLAT edge 
with/without GRIT, (b) FLANGE edge, and (c) ROUND edge 
with/without GRIT. 
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(d) ROUND 



(e) ROUND W/ GRIT 

105 n 



FIGURE 31. Scaled flap edge noise spectra for (a) FLAT edge, 
(b) FLAT edge with GRIT, (c) FLANGE edge, (d) ROUND 
edge, and (e) ROUND edge with GRIT. 


hydrodynamic flow pressure field and the near-field 
noise. The noise generation distribution along the flap 
edge is determined. Two different noise prediction 
methods are developed and successfully validated from 
(1 ) a causality approach, utilizing cross-spectra between 
noise and surface pressures, and (2) an edge-noise 
scatter solution. Both methods use unsteady pressure 
data taken by sensors at the flat flap edge. A CFD flow 
solution is used as a guide for the basic flow pattern 
description, but simplified flow calculations provide the 
necessary boundary layer flow parameter inputs for the 
noise predictions and scaling. The predictions provide 
different but consistent theoretical bases for 
understanding the noise production at the edge. 
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